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We report a complete set of ^^Co NMR data taken on the x — 2/3 phase of sodium cobaltates 
Naa;Co02, for which we have formerly established the in plane Na ordering and its three dimensional 
stacking from a combination of symmetry arguments taken from Na and Co NQR/NMR data. Here 
we resolve all the parameters of the Zeeman and quadrupolar Hamiltonians for all cobalt sites in 
the unit cell and report the temperature dependencies of the NMR shift and spin lattice relaxation 
Ti data for these sites. We confirm that three non-magnetic Co^+ (Col) are in axially symmetric 
positions and that the doped holes are delocalized on the nine complementary magnetic cobalt 
sites (Co2) of the atomic unit cell. The moderately complicated atomic structure resumes then in 
a very simple electronic structure in which the electrons delocalize on the Co2 kagome sublattice 
of the triangular lattice of Co sites. The observation of a single temperature dependence of the 
spin susceptibilities indicates that a single band picture applies, and that the magnetic properties 
are dominated by the static and dynamic electronic properties at the Co2 sites. We evidence that 
they display a strong in plane electronic anisotropy initially unexpected but which accords perfectly 
with an orbital ordering along the kagome sublattice organization. These detailed data should now 
permit realistic calculations of the electronic properties of this compound in order to determine the 
incidence of electronic correlations. 



I. INTRODUCTION 

The cobaltates Na2:Co02 are layered oxyde materi- 
als somewhat similar to the cuprates in as much as the 
charge doping of the C0O2 layers is controlled on a large 
range by variation of the Na content. This can be put 
in parallel with the doping of the cuprates by chemical 
substitutions on the layers separating the Cu02 planes. 
One significant difference with the cuprates is that the 
Co of the C0O2 plane are ordered on a triangular lat- 
tice and not on a square lattice as for the Cu02 plane of 
the cuprates. In this configuration the large crystal field 
on the Co site favors a low spin state in which orbital 
degeneracy influences significantly the electronic proper- 
ties and may yield large thermoelectric effectsii A rich 
variety of other physical properties ranging from ordered 
magnetic states^ high Curie- Weiss magnetism and metal 
insulator transition)^ superconductivity^ etc have then 
been observed on the cobaltates. 

In most cases the theoretical calculations considered 
to explain the physical properties of cuprates and cobal- 
tates have assumed uniform delocalization of the carriers 
in the layers^ and the incidence of the coulomb poten- 
tial of the ionic dopants has been generally considered as 
unimportant. 

Experimentally though, the disorder due to the 
dopants has been shown to have an important incidence 
on the local electronic properties in the cuprates and can 
blur the physical properties and may drive as well sys- 
tem specific, that is non generic properties of the Cu02 
planes. For instance static stripe charge organization of 
the Cu02 planes have been seen only on specific systems 
such as La2-2;Baa;Cu04 for which static distortions lead 
to an orthorombicity of the structure)^ Some anisotropic 
electronic properties of YBa2Cu306-i-a: compounds might 



as well tentatively be associated with the incidence of the 
ordered CuO chains which drive the doping of the Cu02 
planes ii 

In the cobaltates, NMR experiments and structural in- 
vestigations have given evidence that for x ^ 0.5 a large 
interplay occurs between atomic arrangements and elec- 
tronic properties, as the Na are found to be ordered^"— 
It has been shown earlier that this ordering is associ- 
ated quite systematically with cobalt charge dispropor- 
tionation into 3~^ and « S.S"*" Co charge statesJ^ii^"— 
As the cobalt ions in sodium cobaltates are in low spin 
configurations the Co'^'^ has an electronic spin S = 
and appears to be inert magnetically, when compared 
to the other Co sites with higher charge state on which 
holes delocalize (formally Co'*+ should have S = 1/2). 
However the actual Na atomic order, the organization of 
the Co sites and their electronic properties has not been 
so far determined experimentally, except for the case of 
X — 1/2 where such a large charge differentiation does not 
occur Jiii^ This did not therefore permit yet any theoret- 
ical calculations, even in the LDA approximation, based 
on an actual structure. 

We have however established for long that one specific 
structure of Na2:Co02 was quite stable for x « 0.7 and we 
could ascertain recently that it corresponds to a a; = 2/3 
atomic structure which could be determined by combin- 
ing NMR and NQR experiments^^ The Na organization 
in this structure, which agrees with GGA calculations;^ 
consists of two Na on top of Co sites (the Nal sites) and 
six on top of Co triangles (the Na2 sites). This twelve 
Co unit cell and the stacking between planes has been 
confirmed by x-ray diffraction experimentsi^ It results 
in a differentiation of four cobalt sites in the structure, 
two nearly non-magnetic Co^"*" and two more magnetic 
sites constituting a kagome sublattice of the triangular 
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FIG. 1. (Color online) Two dimensional structure of the Co 
planes of Na2/3Co02 compound and the Na environments of 
the four Co sites differentiated in this structure are displayed. 
Cobalts in the Cola and Colb positions have a threefold sym- 
metry and one or two sodium ions in Nal position (red balls) 
as nearest neighbors which drive them into a non-magnetic 
Co"^"*" charge state. In the Co2a and Co2b positions cobalts 
have nearest neighbours sodium ions in Na2 positions (green 
balls). The small difference of local environment of these two 
Co sites only yields negligible distinction between their elec- 
tronic properties. This results in a 2D structure where these 
Co sites form a kagome lattice.— >2i 



Co lattice on which the holes are delocalized (see Fig.[T|). 

Let us point out that a conflicting suggestion concern- 
ing the structure of this phase has been proposed in 
Ref.[22. In Appendix|^we demonstrate that our data are 
quite incompatible with this alternative structure, and 
that nothing allows so far to establish the occurrence of 
oxygen vacancies in our samples. 

In the present paper we shall present a further im- 
portant experimental step which consists in the deter- 
mination of a full set of ^^Co NMR data which, to- 
gether with NQR data, allow us to unravel the electronic 
properties of the four considered Co sites. We confirm 
that the electronic structure resumes in two type of Co 
sites, the non-magnetic Col sites with charge state 3+ 
and the Co2 sites on which holes delocalize. The mag- 
netic properties sensed by all four sites have a single T 
dependent behaviour which confirm a single band pic- 
ture suggested initiallyJ^ We evidence unexpected large 
anisotropics of the magnetic NMR shift and spin lattice 
relaxation of the Co2 sites, which have direct implications 
on the electronic structure of the C0O2 plane. These re- 
sults should trigger realistic calculations to explain the 
correlated electronic structure of these compounds and 
their evolution with hole doping. 

The paper is organized as follows. Sections HHTVI 
contain all experimental aspects of the paper, while in 
Sec. |VI] we collected and discussed all the results which 
are relevant for further theoretical considerations. 

In Sec. [n] we recall briefly the methods used to synthe- 
size single phase samples for x = 2/3 and to orient the 



samples in an applied field. More details on the NMR 
tests of the quality of the sample alignment are given 
in Appendix |B] In Sec. IIIII we recall some basics of the 
NMR spectroscopy of nuclear spins I > 1/2 in presence 
of electric field gradients (EFG) which induce quadrupole 
effects. We report then the ^^Co NMR spectra detected 
in Naa;Co02 samples when the field is applied in both 
directions H \\ c and H J- c, and the experimental proce- 
dures used to separate the contributions of the different 
sites to the NMR spectrum. We describe then the simu- 
lations of the NMR spectra which allowed us to evidence 
the large in plane anisotropy of the NMR shift of the 
magnetic sites. 

In Sec. |TV]we present the T dependencies of the NMR 
shifts of all Co sites in the structure and compare them 
to that observed on the Na NMR (technical details are 
reported in Appendix [C|) . This allowed us to separate 
the T dependent spin contribution to the NMR shifts 
from a T independent term that we assign to an orbital 
contribution. We evidence that both terms are non-axial 
on the Co2 sites. 

In Sec. |V] the technical aspects of the measurements of 
the ^^Co NMR spin lattice relaxation are presented and 
the data obtained on the two types of cobalt sites are ana- 
lyzed by comparing the Co Ti data to that taken formerly 
on ^■^Na. This permits us as well to separate the spin and 
orbital contributions to (TiT)~^. The anisotropics of the 
data for the spin contributions can be explained by those 
obtained for the hyperfine coupling. These results allow 
us as well to demonstrate that the spin and orbital con- 
tributions to (TiT)^^ are at least an order of magnitude 
larger for Co2 than for Col sites. 

In Sec. IVII we then discuss the ensemble of results on 
the EFG, the NMR shifts and {TiT)'^ data. While so 
far the qualitative aspects of the data had lead us to 
consider that the Col sites were Co"^"*", the measured pa- 
rameters allow us to get a little bit further and to obtain 
an upper limit on the charge occupancy of the Col sites. 
More importantly this allows us to consider the origin of 
the anisotropy on the Co2 sites and to conclude that it 
results from the distribution of holes between the axial 
aig orbitals and the in plane e' orbitals on the Co2 site. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL TECHNIQUES 

Experiments on sodium cobaltates Naa;Co02 have 
been performed on both powders or single crystals. Usu- 
ally, NMR experiments, which provide information on the 
bulk of the samples are expected to bring more informa- 
tion on single crystals than in powder samples. However 
this has not be proven to be the case so farj ^^'^'^ Indeed, 
we have shown that homogeneous single phase powder 
samples can be synthesized and can be well controlled 
by x-rays and NMRf^^ Also the penetration of rf mag- 
netic field into conducting cobaltates single crystals is 
limited by the skin effect. Thus for our study we used 



powder samples. However, as they consist of particles 
oriented randomly, the resulting angular distribution of 
quadrupole splittings introduces as well some difficulties 
in the analysis of the NMR spectra. 



Here we should mention that the Stycast perfectly pro- 
tects the powder from environment (water) influencei^ 
For this Na concentration the samples prepared 6 years 
ago did not show any changes in the NMR spectra. 



A. Samples preparation 



III. ^^Co NMR SPECTRA IN THE Naa/sCoOa 



The methods which have been used to synthesize repro- 
ducibly single-phase powder samples with sodium con- 
tent x=0.67 have been reported in Ref. [2ll As described 
there, we used three different routes to synthesize homo- 
geneous single phase samples of Na2/3Co02 compounds: 

(1) direct synthesis from a stoichiometric composition 
of C03O4 and Na2C03. 

(2) from a mixture of cobaltates with calibrated com- 
positions synthesized previously (such as Nai/2Co02 and 
Nao.riCoOa^S) taken in a proper ratio. 

(3) by de-intercalation of Na from Nao.7iCo02 by an- 
nealing it at 700° C - out of its own stability temperature 
range. 

Whatever the synthesis procedure used, the X-ray 
spectra of these Na2/3Co02 samples displayed the same 
c axis parameter and diffraction spectra, including the 
satellite Bragg peaks due to Na order j^ However we no- 
ticed an important difference between these materials re- 
garding the possibility to orient them in a magnetic field. 



B. Powder sample alignment 

As the room T susceptibility of sodium cobaltates is 
known to be anisotropic?^ we have used this anisotropy 
to align the single crystallites of powder samples in a 7 T 
field, by mixing the samples with Stycast 1266 epoxy 
resin which cured in the field. Thus in our samples crys- 
tallites c axes were aligned in the same direction, but ah 
planes of different crystallites are randomly distributed/^ 
To be successful, this procedure requires a powder in 
which individual grains should be single crystallites and 
should not form clusters. Also the shape and packing of 
the grains should not prevent them rotating freely in the 
magnetic field. 

The best aligned samples were obtained from pow- 
ders which were directly synthesized along route (1) with 
grain sizes of 50-100 microns, but a residual unreacted 
C03O4 could not be avoided in many cases and the cor- 
responding ^^Co NMR signal of C03O4 could be seen in 
the spectra. The powders obtained along route (3) also 
usually align well in the magnetic field. Those synthe- 
sized along route (2) were most difficult and sometimes 
impossible to align. In such a synthesis the small grains 
are probably randomly welded to each other, so that the 
final powder does not consist of single crystallites. 

For a perfect alignment the samples should look like 
a single crystal in the c direction but with full ah plane 
disorder. Details on the NMR tests of the quality of the 
sample alignment are given in Appendix [Bl 



A. NMR Background 

Generally in solid state NMR an atomic nucleus with 
spin / and quadrupole moment Q has its spin energy 
levels determined by the Zceman interaction with a ex- 
ternal magnetic field Hq and the quadrupolar interaction 
with the electric field gradient (EFG) on the nucleus site. 
Therefore the Hamiltonian consists of two parts - Zeeman 
Hamiltonian "Hz and quadrupolar Hamiltonian "Hq, and 
can be written22i2^: 



H = Hz+Hq^ -7/1/(1 + K)Ho + 



eQ 



2I{I- 1 



-IVI, (1) 



where 7 is gyromagnetic ratio. 

The physical properties of the studied compound are 
hidden in the two tensors in this Hamiltonian - magnetic 
shift tensor K and the EFG tensor V . The principal axes 
of both tensors are associated with the local structure and 
assuming that the principal axes {X, Y,Z) of both tensors 
coincide, the Zeeman and quadrupolar Hamiltonians can 
be re-written: 



T-Lz = -7/1^/0(1 + Kaa)Hoa; 



-Hn^ 



eQ 



l\ Z^ ^aala 



(2) 



21(1 

-^[34-/(/-l)]+K-.y)(/i-4)), 

where a ~ X,Y, Z, and the three quadrupolar frequen- 
cies 



" [2/(2/ - l)h] 

are linked by Laplace equation J^a ^aa = X^a ^a = 0. 
Therefore it is more common to use two parameters - the 
quadrupolar frequency vq = vz corresponding to the 
largest principal axis component Vzz of the EFG tensor 
and the asymmetry parameter t\ = {Vxx — Vyy)/Vzz 
(here the principal axes of the EFG tensor are chosen 
following \Vzz\>\Vyy\> \Vxx\)- 

The nuclear spin of ^^Co is / = 7/2 and therefore for 
a given direction of applied magnetic field Hq relative 
to the crystallite the NMR spectrum for a single Co site 
consists of 7 lines - a central line which corresponds to the 
— 2 "!->■ 2 transition and 6 satellites corresponding to the 
other 771 f-> (m — 1) transitions. The position of the cen- 
tral line is determined by the applied field and the values 
of the magnetic shift K and second order quadrupolar 
perturbationj^ The distance between the two satellite 



lines mo (m — 1) and — m o — (m — 1) depends on 
the orientation of the external field with respect to the 
principal axis of EFG tensor described by the spherical 
angular coordinates 9 and f and can be expressed asi^ 



Ai^ = Vq{7 



1 



)(3cos2 6l- l + ?7sin26icos2(^). (3) 



B. Experimental techniques 

The NMR measurements were done using a home-built 
coherent pulsed NMR spectrometer. NMR spectra were 
taken "point by point "with a 7r/2 — r — 7r/2 radio fre- 
quency (rf ) pulse sequence by varying the magnetic field 
in equal steps. The minimum practical r values used in 
our experiments was 7 ^s. The usual 7r/2 pulse length 
was 2 /is. 

Experimentally the value of the magnetic shift (in %) 
for a given point in the spectrum either in frequency do- 
main or in the field domain can be calculated using: 



K = 



'■'ref 



Vref 



Bref — B 

B ' 



(4) 



where Vref and Bref are reference frequency and field 
values, respectively, which are connected by the relation 
Vref = {l/'^T^)Bref- In this work ^•^Na shifts are given 
with respect to the ^'^Na resonance in a NaCl water solu- 
tion. The reference for '^^Co shift is given for the 7/27r — 
10.054 MHz/T value. 



C. Decomposition of the spectra in slow and fast 
relaxing components 

The raw ^^Co NMR spectra in the Na2/3Co02 com- 
pound are complicated - as an example we show in Fig. [5] 
the NMR spectra measured in a well oriented sample at 
T =5 K for two directions of the applied magnetic field 
Hi). Former ^^Co NMR dataM taken on this phase have 
already allowed us to evidence distinct Co NMR lines, 
but those NMR spectra were somewhat difficult to ana- 
lyze fully, as one needs to determine altogether the EFG 
parameters and NMR shifts of the various sites. The 
experiments were furthermore complicated by the need 
of a quasi perfect alignment of the powder sample with 
respect to the applied magnetic field. 

We foundi^i2i that some Co sites have both very short 
nuclear spin-spin (T2) and spin-lattice (Ti) relaxation 
times - we denote them as Co2 type. The Co sites which 
have much longer relaxation times could be also isolated 
- we call these cobalts as Col type. 

In Fig. [5] we show a decomposition of the experimental 
spectra in the Col and Co2 contributions. The spectrum 
measured with r = 7/is between rf pulses contains con- 
tributions of the NMR signals of both type (Col-|-Co2) 
of ^^Co nuclei in the sample, whereas for r — 200/is the 
only remaining contribution to the signal is that from 
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FIG. 2. (Color online) Examples of ''''Co NMR spectra of an 
oriented sample of the Na2/3Co02 compound taken with the 
applied magnetic Ho _L c (upper panel) and Ho \\ c (lower 
panel). In each panel the top spectra are taken using a short 
time interval r = 7/is between rf pulses. They involve the 
NMR signals of all cobalts in the sample. The signals of the 
slow relaxing Col sites and ^^Na are isolated in the second 
spectra taken with r — 200/is. The NMR spectra of the fast 
relaxing Co2 sites could be obtained either by subtraction of 
short-r and long-r spectra (Co2 by difference), or by satu- 
ration of the slow relaxing Col and ^^Na NMR signals by 
an additional rf pulse (Co2 by saturation). Although some 
remnants of the Na and Col signals remain in the difference 
spectra, the actual quadrupole splitting of the Co2 spectrum 
is quite well resolved for Ho \\ c. 



the slow relaxing ^^Col nuclei. One can then obtain the 
spectrum of the fast-relaxing part Co2 - by subtracting 
the rescaled slow-relaxing part of the spectra from the 
short-T spectra. The rescaling can be done empirically 
by trial and error, or using an estimation of the T2 decay 
for the slow relaxing spectrum. In the Fig. [3] the spec- 
trum obtained so for the fast-relaxing part Co2 is shown 
as " Co2 by difference" . 

Slow and fast-relaxing cobalts can be separated as well 
by the large difference in their spin-lattice relaxation 
timeS)^ It is possible to suppress (or more correctly, to 
saturate) the signals of ^'^Na and slow-relaxing Col type 
sites by using an additional 7r/2 pulse with some delay 



(usually 400-1000 fis) before the 7r/2 - r - 7r/2 pulse 
sequence. In the spectrum obtained after such a pulse 
sequence the intensities of the slow-relaxing cobalt lines 
is reduced and the lines of the fast-relaxing cobalts are 
better resolved. Spectra of Co2 obtained by this method 
are shown in Fig. [2] as " Co2 by saturation" . 

The shapes of the fast-relaxing Co2 NMR spectra ob- 
tained by the two methods are not exactly the same as 
one can see in the Fig. [21 The compliance between the 
two methods is much better when the external magnetic 
field is applied parallel to the c axis of the sample. The 
main reason for the discrepancy is that the nuclear re- 
laxation rates differ in the different parts of the NMR 
spectra. Therefore it is not so easy to subtract fully or to 
suppress fully the Col spectra, so that some remanent of 
the intense Col lines remains in the "Co2 by difference" 
spectra. Also the overlap with the ^■^Na NMR signal 
creates additional difficulties - one can see the residual 
traces of the sodium signal in the obtained Co2 spectra. 
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D. Slow relaxing Col NMR spectra 

Figures [3] and U show spectra of the slow relaxing 
Col at two temperatures, 5 K and 60 K, respectively. 
These experimental spectra are quite simple as in both 
directions of applied magnetic field Hq J- c and Hq \\ c 
they are nearly symmetric with respect to the center of 
the spectra. The largest splitting between satellites for 
slow relaxing Col is observed when the applied magnetic 
field Hq is parallel to the c axis of the sample. This fact 
demonstrates that the main principal axis Z of the EFG 
tensor for Col in the Na2/3Co02 compound is parallel 
to the crystallographic c axis. The fact that the distance 
between outer singularities in both directions Hq _L c and 
Hq II c scales as 1:2 shows that these cobalt sites have an 
axial EFG tensor with asymmetry parameter 77 close to 
zero (see Eq. [S]). 

The ^^Co NMR spectra of the slow relaxing Col can 
be described by only 2 unequivalent cobalt sites which 
we assign to the Cola and Colb sites of Fig. [T] In figures 
[3]and [3]we show computer simulations of the ^^Co NMR 
spectra that wc could perform. As one can see the agree- 
ment between the simulated spectra and the data is very 
good. These simulations of the NMR spectra allow us to 
deduce vq , 77 and the magnetic shift Kaa of the central 
line for the two Col sites, which are summarized in the 
Table D 

These results are in good agreement with the parame- 
ters deduced from ^^Co NQR spectra which also display 
two slow relaxing sites.— Moreover, the simulation of the 
^^Co NMR spectra also confirm the ratio 2:1 between the 
lines intensities corresponding to the Cola and Colb sites 
as can be seen in Fig. |3]and |H 



FIG. 3. (Color online) '^''Co NMR spectra ((red) dots joined 
by a full line) of the slow-relaxing Col sites measured at 
T =5 K in the two directions of the applied magnetic field 
Ho which were shown in Fig. [2l The simulation of the spec- 
tra of Cola site are shown by (magenta) dash-dotted lines, 
those of the Colb spectra are shown by (cyan) dashed lines. 
The solid black line corresponds to the sum of the Cola and 
Colb simulated spectra. 



E. Details of the simulations of the NMR spectra 

The simulation procedure used to reproduce the pow- 
der NMR spectra requires an averaging of computed sig- 
nals for all possible orientations of the powder particles. 
In our samples the c axes of the crystallites are aligned 
but the a or 6 axes are at random. To simulate a pos- 
sible imperfect alignment of the powder particles, we in- 
troduced in our simulation a distribution of crystallite 
c axes orientations which we described by a Lorentzian 
function. For well oriented samples (like in Fig. [2]) the 
deviation of c axis orientations was within ±5°. Such an 
angular distribution described very well the ^'^Na NMR 
spectra (see Appendix [B)) and has been used then in the 
simulations taken for all cobalt sites in a given sample 
(see figures [Ml]). 

Simulations of the spectra were done taking into ac- 
count not only the splittings of the energy levels de- 
scribed by the Hamiltonian of Eq. [2] but also the tran- 
sition probabilities due to the radio frequency field exci- 
tation using the algorithm described in Ref. [2^. Also to 
simulate the broadening of the NMR lines we used a tri- 
angular filtering function with bandwidth 0.15 MHz for 
both Col and Co2 sites at T=60 K. At T=5 K the larger 
broadening of the Co2 site spectrum required an increase 




TABLE I. Parameters used in Hamiltonian [2] for com- 
puter simulations of the ^^Co NMR spectra shown in fig- 
ures |3]|6] for different cobalt sites in Na2/3Co02. The prin- 
cipal components of the tensors are given in the princi- 
pal axes {X, Y, Z) of the EFG tensor. Site occupancy 
ratios Cola/Colb/Co2a/Co2b=l/2/3/6 were used in the 
simulations.— 
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FIG. 4. (Color onhne) Experimental ^^Co NMR spectra 
((red) dots joined by a full line) of the slow-relaxing Col sites 
measured at T =60 K in the two directions of the applied 
magnetic field Hq. The simulation of the spectra of Cola site 
are shown by (magenta) dash-dotted lines, those of the Colb 
spectra also shown by (cyan) dashed lines. The solid black 
line corresponds to the sum of the Cola and Colb simulated 
spectra. 



of this bandwidth to 0.5 MHz, while 0.15 MHz could be 
kept for the Col sites. 



F. Fast relaxing Co2 NMR spectra 

Figures [5] and [5] show the spectra of the fast relaxing 
Co2 at two temperatures, 5 K and 60 K, respectively. De- 
spite the difficulties in obtaining fast relaxing Co2 NMR 
spectra, some facts can be settled. First of ah the largest 
splitting between satellites for Co2 is observed for Ho \\ c. 
Therefore for Co2 the main principal axis Z of the EFG 
tensor is parallel to the crystallographic c axis in the 
Na2/3Co02 compound. 

The Co2 spectrum for iJo || c at T =60 K (see Fig. ^ 
is quite simple, well resolved and clearly evidences two 
inequivalent Co2 sites with different values of quadrupolc 
frequency vq - we denote these sites as Co2a and Co2b. 
The relative intensities of the NMR signals of these sites 
corresponds to a ratio Co2a:Co2b close to 1:2. These 
facts are in very good agreement with the ^^Co NQR 
spectra which also display two fast relaxing sites with a 
similar intensity ratio i^ 

At low temperature T =5 K the Co2 spectrum in the 
Hq 11 c case (see Fig. [5]) is less resolved due to the large 







Cola 


Colb 


Co2a 


Co2b 




Kxx 


1.93 


2.08 


9.52 


10.4 




Kyy 


1.93 


2.08 


4.36 


4.76 


NMR 


Kzz 


3.48 


3.20 


3.54 


3.80 




^ol 


1.19 


1.386 


2.276 


2.591 


5K 


V 








0.362 


0.358 




uxl 


-0.595 


-0.693 


-0.73 


-0.83 




vy:_ 


-0.595 


-0.693 


-1.55 


-1.76 




Kxx 


1.98 


1.98 


5.19 


6.0 


NMR 


Kyy 


1.98 


1.98 


2.80 


3.31 




Kzz 


2.67 


2.72 


2.67 


3.72 


60 K 


^^ 


1.165 


1.340 


2.19 


2.55 




V 








0.362 


0.358 


NQR 


ucA 


1.193(1) 


1.392(1) 


2.187(1) 


2.541(1) 


4.2 K 


V^ 


<0.017 


<0.016 


0.362(5) 


0.358(4) 



^ Since in the NMR the sign of the EFG cannot be determined, 
the sign of values of Ua could be reversed. 



magnetic broadening of the lines, but the satellites of the 
two sites could still be seen. 

If we compare now the Col and Co2 NMR spectra 
measured in the Hq -L c it is easy to see that Co2 spec- 
tra are asymmetric with a larger intensity in the lower 
fields. At T=60 K some lines are resolved with clear sin- 
gularities but at T =5 K it becomes nearly impossible to 
resolve the different lines and we could only detect very 
weak singularities. As the orientation procedure of our 
samples is such that the ab planes of different crystallites 
are randomly distributed then the experimental spectra 
in the Hq _L c direction are powder spectra. This ex- 
plains why they are less resolved than the Hq \\ c spectra. 
However some singularities could still be perceived even 
in the broadened spectrum at T =5 K (Fig. [5]). Such a 
complexity is apparent as well in the data taken by many 
others groups in sodium cobaltates with sodium content 
X sa 7, 30,31 Comparing our oriented powder NMR spec- 
tra with NMR data taken on single crystalsi^i^ allows us 
to conclude that twinning of the ab plane is quite com- 
mon in those single crystals. 

Therefore the analysis of a single Hq _L c spectrum of 
Co2 appears quite difficult. Independently of the method 
used to isolate the fast relaxing Co2 spectra we do how- 
ever find that the positions of the different singularities 
in the Co2 spectra are the same as one can see in Fig. [5] 
As usual such singularities in a powder NMR spectrum 
appear at positions where the two principal axes X and 
Y of the EFG and magnetic shift tensor in the ab plane 
correspond to extremal positions of the NMR signal for 
randomly distributed crystallites. Furthermore, as was 
shown in Ref. [ij the sodium and cobalt sites in the 
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FIG. 5. (Color online) ^®Co NMR spectra ((green) dots joined 
by a full line) of the fast-relaxing Co2 sites measured at 
T =5 K in the two directions of the applied magnetic field 
Hq. Those were obtained by saturation of slow-relaxing Col 
and Na NMR signals and are taken from Fig. [J] The sim- 
ulation of the spectra of the Co2a site are shown by (red) 
dash-dotted lines, those of the Co2b spectra also shown by 
(blue) dashed lines. The solid black line corresponds to the 
sum of the Co2a and Co2b simulated spectra. The large devi- 
ations seen for B~6.45 Tesla and B«5.9 Tesla are due to the 
imperfect saturation of the Col and Na contributions. 



FIG. 6. (Color online) '^^Co NMR spectra ((magenta) dots 
joined by a full line) of the fast-relaxing Co2 sites measured 
at T =60 K in the two directions of the applied magnetic field 
Ho . Those were obtained by subtraction of long-r spectra 
from the short-r spectra. The simulation of the spectra of the 
Co2a site are shown by (red) dash-dotted lines, those of the 
Co2b spectra also shown by (blue) dashed lines. The solid 
black line corresponds to the sum of the Co2a and Co2b sim- 
ulated spectra. The large peaks between 6.2 and 6.4 Tesla are 
remnants of the Col NMR signals, incompletely subtracted. 



Na2/3Co02 compound feels a common T dependence of 
the spin susceptibility, which we could determine directly 
by measuring the ^'^Na magnetic shift ^^K{T)^ We did 
then study the evolution of the positions of the singular- 
ities in the Co2 spectrum versus temperature. By plot- 
ting these singularities versus the spin susceptibility as 
monitored by ^^K ^ the temperature being then an im- 
plicit parameter, we could sort out those singularities 
which correspond to Co2 central lines and those which 
can be assigned to quadrupolar satellites, as detailed in 
Appendix O 

This allowed us to resolve then fully the EFG and shift 
tensors of Co2a and Co2b sites and to understand that 
not only the EFG tensor but also the shift tensors arc 
highly asymmetric. We could then succeed to perform a 
full simulation of the '^^Co NMR spectra of the fast relax- 
ing cobalts in the two field directions, as shown in figures 
[S] and m with the parameters collected in the Table ID 
There we also reported the parameters measured for the 
Cola and Colb sites, as obtained from the simulations 
done in Fig. [3] and Fig. 2] As one can see the obtained 
values of the quadrupolc frequency vq and of the asym- 



metry parameter 77 for all sites are in good agreement 
with the values obtained by NQR<^ 



IV. T DEPENDENCE OF THE COBALT 
MAGNETIC SHIFTS 

The principal components of the magnetic shift tensor 
for the cobalt site i in a direction a consist on three parts 
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Here iCf^ is the chemical shift due to the diamagnetic 
susceptibility of the inner shell electrons, K™^ reflects 
the on site orbital effect of the valence electrons, propor- 
tional to the orbital part of the electronic susceptibility 
Xi^a which is usually temperature independent. The last 
term, the spin shift K^^^^, is proportional to the local 
electronic spin susceptibility xt aC^)- 

The remarkable feature of the Na2/3Co02 phase of the 
sodium cobaltates is the large low T variation of this 



last term, which, as already pointed out above, allowed 
us to resolve the Co NMR spectral The ^^Na NMR 
has negligible chemical and orbital NMR shifts and is 
largely dominated by the T dependent spin term which 
has been found nearly isotropicjSi So the comparison of 
^^Ki^a with the isotropic contribution ^^Kiso also allowed 
us to separate the orbital and spin contributions to the 
cobalt magnetic shift, as detailed below. 



Col shifts 



TABLE II. ^^Co magnetic shift parameters obtained for the 
4 cobalt sites. 



Site 






KTx(%) 



Cola 



Colb 



-0.2(2) 
-0.2(2) 
4.0(2) 



0.2(2) 
0.2(2) 
3.5(2) 



1.92(2) 
1.92(2) 
1.86(5) 



1.92(2) 
1.92(2) 
1.86(5) 



Co2a 



17.9(4) 
7.1(3) 
4.6(2) 



3.1(2) 
2.1(2) 
1.92(5) 



Co2b 



19.9(4) 
7.9(3) 
5.2(2) 



3.3(2) 

2.3(2) 
1.95(4) 



In Fig. [7^ we show the T dependence of the magnetic 
shifts of Cola and Colb sites which arc plotted in Fig.[7lD 
versus the NMR shift of ^-^Na. For Cola and Colb sites 
the Z-components of the magnetic shift arc quite similar, 
with a substantial T dependence which is proportional to 
^^K at temperatures below 150 K. 

The in-plane X and Y components of the shift are 
indistinguishable for Cola and Colb sites, and their shift 
value does not change from 150 K down to T =40 K. 
Below this temperature some small difference in the in- 
plane shifts for Cola and Colb sites appears. This was 
also seen by others;^ which points out that this effect 
is somewhat characteristic of this x = 2/3 phase of Na 
cobaltates. But as one can see in Fig. [71 these distinct 
variations of NMR shifts of Cola and Colb are quite 
small when compared with the T variations of the NMR 
shifts of both Col sites for Hq || c. So we can consider so 
far that the in plane shift of the Col sites is practically 
T independent below 150 K. 

Using Fig. [7}3 one can get the slopes of the linear fits of 
^^K versus ^^K which give the relative magnitudes of the 
spin contributions K^ ^f^^K, while the ^^K=0 intercepts 
give the estimates of the T independent contribution to 
the ^^Co NMR shift. The values found here, reported 
in Table [ill arc much larger than those expected for the 
chemical shift Kf^^ < 0.1%. They are obviously domi- 
nated here by the Co orbital NMR shift K°^^, so that we 
have neglected Kf'^^ in all our analyzes. 

But, as seen in Fig. [71 significant increases of ^^K^ 
with respect to the low T linear ^^K dependence are ob- 
served for T > 150 K, that is for ^^K < 0.08. Those in- 
creases are associated with the onset of Na motion which 
has been detected at ^200 K from ^'^Na data in Ref. QQ. 
There we have discussed that this Na motion induces a 
significant Col-Co2 site exchange which induces such an 
increase of Col NMR shifts, while the charge dispropor- 
tionation already occurs above room temperature. 



B. Co2 shifts 

In Fig. [8^ we show the T dependencies of the mag- 
netic shifts of Co2a and Co2b sites and in Fig.[BjD they are 
shown versus ^^K. In this figure we put some experimen- 
tal points which were determined accurately at specific 



temperatures, but due to the relatively small number of 
Co2a sites the intensity of its NMR signal is weak. For 
this site we plot lines corresponding to the Kx and Ky 
components of its shift obtained from the analysis de- 
tailed in Appendix B. 

As one can see in Fig. [BJd a linear dependence is found 
for all components of the magnetic shift for both Co2a 
and Co2b sites below 120 K in a range where no sign of 
Na motion could be seen. So this figure in combination 
with Fig. [TJd confirm that all 4 cobalt sites pertain indeed 
to a unique sodium cobaltate phase in which a single T 
variation characterizes the local xtaC^)- ^^ collected 
in Table [ll] the relative magnitude of the spin contribu- 
tions K^ ^f^^K and K""^^ for all three components of the 
magnetic shift for the Co2a and Co2b sites together with 
those of the Col sites. 

While the magnetic shift tensors are axial for the Cola 
and Colb sites, we see in Table [IT] that both Co2a and 
Co2b sites demonstrate strong in-plane anisotropy of 
both the orbital shift and the spin hyperfine coupling. 
This totally unexpected anisotropy is a quite important 
observation achieved in the present investigation. In our 
previous paper concerning the ^^ Co NMR in this phase^^ 
as well as in all the existing publication a^^i^^'^°i'^^ such 
anisotropy was not anticipated. This had consequently 
led us to invoke a third group of axial cobalt sites Co3 
with large shift values. In appendix [B] we detail why im- 
perfect alignment of the sample powder had led us as well 
to overlook initially this anisotropy. 



V. ^^Co SPIN LATTICE RELAXATION 

When thermal equilibrium of the nuclear spins is dis- 
turbed by rf pulses, the equilibrium nuclear magnetiza- 
tion is recovered by various relaxation processes which re- 
fiect the interactions in the spin system and the magnetic 
and electronic properties of the materials. The recovery 
of the nuclear magnetization along the applied field oper- 
ates through nuclear spin lattice relaxation (NSLR) pro- 
cesses characterized by the spin-lattice relaxation time 
Ti or by the nuclear spin- lattice relaxation rate 1/Ti. 

In systems with unpaired spins the dominant Ti pro- 
cess is due to local field fiuctuations induced by the dy- 
namics of the local electronic magnetization. Theoreti- 
cally, the spin contributions to {TiT)~^ may be written 
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FIG. 7. (Color online) (a) T-dependencies of the ^^Co NMR 
shifts of the Cola (filled (magenta) symbols) and Colb (open 
(cyan) symbols) sites for two directions of the external mag- 
netic field Ho- For Ho \\ Z (circles) the Cola and Colb shifts 
exhibit the same T dependence with a small difference in ab- 
solute values. In the Ho J- Z direction for both sites Cola 
and Colb the T dependence is much weaker and a small dif- 
ference between the shifts of these sites appears only below 
T—40 K. For both Cola and Colb sites the magnetic shifts 
are axial, (b) The data of (a) are plotted versus the ^"^Na 
NMR shift showing linear slopes at low temperatures and de- 
viations from this linear dependence above 100 K. Such low-T 
linearity allowed us to determine the orbital Kf^^ and spin K^^ 
contributions to ^^K for all three directions i = X,Y, Z - see 
Table HI] 
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FIG. 8. (Color online) (a) T-dependencies of the ^^Co NMR 
shifts of the Co2a ((red) closed symbols and (red) dashed 
lines) and Co2b ((blue) open symbols and (blue) dashed lines) 
sites for the three directions of the external magnetic field Ho 
parallel to the X, Y, Z principal axes of the EFG tensor. The 
method used to determine the data for Ho \\ X,Y directions 
are detailed in Appendix [C] (b) The data of (a) are plotted 
versus the ^"^Na NMR shift. The linear variation of ^^K versus 
^'^ K allowed us to determine the orbital K°i^ and spin Kl^ 
contributions to the ^^''^ K for the three directions i = X,Y,Z 
- see Table HH 



using the imaginary part of the Co dynamical electron 
spin-susceptibility x"(q, fn) as 



(TiT)- 



2%'.fc 
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l^(q)l 



.X"(q, '^n) 
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where A(q) = Aq -f ^A^e"'''' is the wave- vector q de- 
pendent hyperfinc form factor which depends on the hy- 
perfine interaction between the observed Co nuclear and 
electron spins at the same site Aq and the hyperfine in- 
teraction with electron spins at nearest neighbor Co sites 
A M 

Our measurements of the ^■^Na NSLR temperature de- 
pendence allowed us to establish that in sodium cobal- 



tates with sodium content x > 0.6 the dominant correla- 
tions in the electron spin system are ferromagnetic, while 
for lower sodium contents antiferromagnetic correlations 
take overiiSii^ 

Equation [5] applied to different nuclei tells us that 
NSLR data should be similar for all nuclei in the ab- 
sence of any cancelation at some specific q vector due to 
the form factor associated with the geometrical position 
of the nuclear site with respect to the magnetic atoms. 
So the comparison of ^^Co and ^^Na NSLR could help us 
to better characterize the dynamic susceptibility. We re- 
port then here the results of a detailed study of the ^^Co 
NSLR in this Na2/3Co02 phase. First we demonstrate 
the consistency and validity of our measurements, and 
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then compare the results for the different Co sites. We 
follow as well the procedure used for the shift measure- 
ments and compare directly the Co and Na Ti data. 



A. Experimental techniques 

To study the nuclear spin-lattice relaxation process of 
the ^®Co we have used the usual magnetization inversion 
recovery method with three pulses: it — t — 7r/2 — r — tt. 
In this sequence the first pulse rotates the magnetization 
by 180° and the longitudinal magnetization recovered af- 
ter time t is measured with a spin echo sequence with 
interval r. The dependence of the spin-echo intensity on 
delay time t allows to monitor the recovery of the nuclear 
magnetization associated with a given NMR transition: 



M{t)^Ma{l-B-R{t)). 



(7) 



Here Mq is the thermal equilibrium value of magnetiza- 
tion and the parameter B characterizes the actual mag- 
netization after the first pulse at i = (the imperfection 
of the experimental conditions gave typical values B ~ 
1.8 rather than B = 2 expected for a perfect tt pulse). 

The shape of the relaxation function R{t) depends 
on the nuclear transition sampled. For a two-level nu- 
clear system (like the / = 3/2 NQR case) this process 
is exponential R{t) — exp (— i/Ti) and allows a sim- 
ple experimental determination of Ti.— But generally 
for / > 1/2 the nuclear energy levels are differentiated 
by the quadrupole interaction with the crystalline elec- 
tric field (Eq. [1]). For very broad NMR spectra as those 
considered here, the applied first pulse only inverts the 
populations of some levels which, for a given frequency, 
can be selected by properly choosing the applied field. 
Consequently the difference in population between ad- 
jacent levels which are probed by the rf pulses depends 
on the populations of the levels which are not hit by the 
rf pulses. Therefore the magnetization recovery becomes 
multi-exponential 



R{t) =Y^a,exp(-jr] 



(8) 



but is still characterized by a single Ti value. 

The spin-lattice relaxation could be driven either by 
magnetic or quadrupolar fluctuations which correspond 
to distinct transition probabilities and A^ values. How- 
ever in the sodium cobaltates the magnetic relaxation 
mechanism dominates at least at low T^ Using Ref. [Sj 
and|3Jwe have calculated for our values of 77 the Oj and A^ 
parameters for the theoretical relaxation functions R{t) 
for the case of magnetic relaxation by weak fluctuating 
magnetic fields for the different ^^Co transitions - see Ta- 
ble Illl 



TABLE III. Parameters \i and ai of the theoretical spin- 
lattice relaxation function (|S} for different transitions of ^'"'Co 
for the case of magnetic relaxation by weak fluctuating mag- 
netic fields (based on Ref. [33 and[3J). 



Ai 28 

Transition ai 



21 

a2 



15 

13 



10 

04 



6 

a 5 



3 



1 

a? 



^yiT+yimi 0.206 0.068 0.012 

±1 ^ ±1 0.457 0.371 0.001 0.117 0.030 0.012 0.012 

±1 o ±1 0.114 0.371 0.366 0.081 0.008 0.048 0.012 

±1 o ±1 0.009 0.068 0.206 0.325 0.273 0.107 0.012 



B. '' Co spin lattice relaxation results 

We have already seen by NQR that Ti data are iden- 
tical for the Cola and Colb sites. They are identical as 
well for Co2a and Co2b sites but much shorter than on 
the Col sites, which allowed us to confirm that the Co2 
sites were indeed the magnetic sites.— 

Due to the close values of the magnetic shift for all 
4 cobalt sites in the Z direction (see Table HI]) . to avoid 
cross-relaxation between sites it is of course better to take 
data on quadrupolar satellite transitions of different sites 
which are quite separated in frequency (or field). So we 
have measured the spin-lattice relaxation curves on the 
outer satellites -|-| O- -|-| for the sites with larger signal 
intensity, that is Cola for the slow relaxing sites and 
Co2b for the fast relaxing sites. In Fig.[9l we display such 
examples of experimental spin-lattice relaxation curves 
for Cola and Co2b sites. One can see there that they can 
be fitted quite well by the function [8] with the coefficients 
from Table IIIII It can also be seen in Figure [9^ that the 
same Ti value is deduced from the data for the f <-^ | and 
I O I transitions for the Cola nucleus, demonstrating 
convincingly that the data analysis is reliable and that 
the relaxation is indeed magnetic. 

We could also evidence in Fig. [9]3 that the Ti data 
depend of the orientation of the applied field for Co2, 
contrary to the isotropic results found for the Col sites 
in Fig. Hi. 

Similarly, in Fig. llOh we plotted then the T variation 
of the Co2b T^^ data measured in two distinct applied 
fields, showing that Ti values measured are identical, 
which confirms the reliability of the measurements, in 
view of the large spectral changes which occur with in- 
creasing field. We also compare there the data taken on 
Cola with those on Co2b. The Ti values differ by a 
factor as large as 20 which confirms the strong differen- 
tiation already found in previous observations done by 
NQR at low r3i 

Generally, as given in Eq. [S] the electron spin contribu- 
tion to the spin-lattice relaxation rate T^ is governed 
by the fluctuations of the transverse components of the 
effective magnetic field induced at the nucleus by the 
electronic magnetization. In NMR the quantization axis 
is given by the external magnetic field and for the ^^Co 
Tj~ the dominant term is the fluctuating hyperfine flcld 
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FIG. 9. (Color online) Spin-lattice relaxation curves measured 
at T = 80 K for two orientations of the magnetic field, (a) for 
Cola sites the data shown have been taken on two satellite 
transitions | •<->■ | and | ^-^ f . The fits of the data with 
the magnetization relaxation functions ([SI, with parameters 
from the Table [Till give a unique Ti value, independent of the 
field orientation, (b) For Co2b similar good fits are obtained 
for the data, but the relaxation time depends on the field 
orientation. 



whose magnitude is proportional to the hyperfine field 
|A(q)| and local spin susceptibility x"(Q: ^n)- The NMR 
shift and Ti data taken on ^^Na have given evidence that 
these electronic susceptibilities are nearly isotropic r'li^ so 
the spin contributions to 1/Ti on all nuclei are expected 
to scale with each other and the Ti anisotropies should 
be governed by those of the hyperfine couplings |^(q)|. 
For instance 
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Therefore as the spin component of the magnetic shift 
Kxx for Co2 is much larger than that for Col (see Ta- 
ble |TT| it is quite natural to find much shorter Ti val- 
ues for Co2 than for Col. Also, as the components 
of the magnetic shift tensors for Co2 sites are highly 
anisotropic, with a dominant value for Hq \\ X, the trans- 
verse fluctuating local fields are expected much larger 
for Tj^ than for T^^ in qualitative agreement with the 
anisotropy of the Ti data displayed in Fig. ^p. 

However, for Col a distinct situation occurs. As one 
can see in Fig. [9^, for Cola the spin-lattice relaxation is 
isotropic while the spin component to the shift is strongly 
anisotropic - Kzz is much larger than the in plane com- 
ponents which are vanishingly small - see Fig. [7] and Ta- 
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FIG. 10. (Color online) (a) Temperature dependence of the 
NSLR rate T^^ measured in _ffo || c orientation for the Cola 
(magenta empty circles) and Co2b sites (blue closed circles 
and squares, taken for two different applied fields), (b) The 
(TiT)"^ of ^^Co of (a) are plotted versus {T^TY^ values of 
^■^Na, the temperature T being then an implicit parameter in 
the plot. For both Cola and Co2b sites the variations are 
found linear, which allows us to separate the spin and orbital 
contributions to the ^^Co (rir)"^(see text). 



bicini 

To explain that, we are lead to consider whether the 
spin components of the electronic dynamic susceptibility 
do dominate the relaxation, or is there another indepen- 
dent contribution to l/Ti, such as for instance an orbital 
contribution to the spin lattice relaxation? Indeed, usu- 
ally for transition metals two additive contributions to 
T\ are expected 
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So we tried to separate such components similarly to 
what we had done before for the cobalt shifts - see 
Section |lVl Indeed, as was shown in Ref. [13 in the 
Na2/3Co02 compound the ^'^Na spin-lattice relaxation 
is largely dominated by the T dependent spin term while 
any orbital contribution is negligible on that nuclear site. 
So we can plot the values of (TiT)"^ of ^^Co versus 
(TiT)"-^ values of ^-^Na, as done in Fig.fTUb. Good linear 
relations are found for both sites, that is 



'(TiT), 



23 



(TiTj-^ 3.8(1) + 2.4(1) 
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for Cola and 



59/ 



{TiT)z = ^-^(TiT)-!* 270(3) + 23(2) 

for Co2b. This analysis allows indeed to separate in the 
Ti data the spin part from an orbital TiT = constant 
parti^ Let us notice that the ratio of the spin contribu- 
tions for Co2 and Col is about 270/3.8 w 70. Such a 
large ratio compatible with that obtained by NQR— is 
quite expected as the spin part of the Col shift in the X 
and Y directions is practically zero - see Table [ITl 

As for the orbital contribution, one can see in Fig. [TU] 
that it is small with respect to the spin contribution for 
the Co2 site while it becomes dominant for the Col site. 
For instance at T=80 K the Ti is totally dominated by 
the orbital term for the Col site, which perfectly justifies 
the absence of detected Ti anisotropy on the Col site 
at this temperature. Let us point out however that the 
orbital term is still one order of magnitude larger for Co2 
than for Col. 

Overall these Ti data are totally compatible with the 
anisotropy of the Co2 shifts reported in Sec. IIVBI Wc 
shall discuss further the implication of these results on 
the electronic properties of this Na2/3Co02 phase in the 
following discussion section. 



VI. DISCUSSION 



with the imperfect alignment of the powder grains for 
some samples, which had initially lead us to conclude 
that the two other Co sites displayed quite distinct nearly 
axial NMR shiftsJ^ The analysis of the extensive set of 
^^Co NMR spectra done for the two fast relaxing sites 
Co2a and Co2b allow us to establish that these sites have 
nearly identical spin shifts tensors which are unexpect- 
edly totally non-axial, the larger NMR spin shifts being 
in a privileged X direction in plane. The large anisotropy 
of the spin contribution to the 1/Ti data agrees perfectly 
with this NMR shift anisotropy, as Tiz is found shorter 
than Tix- 

We therefore established here that these fast relaxing 
Co2a and Co2b sites have similar anisotropic orbital and 
spin shifts and almost identical Ti valuesr^^ confirming 
that they contribute identically to the electronic prop- 
erties and are dominantly responsible for the magnetic 
properties of this compound. 

This allows us then to conclude that although the lo- 
cal environment of Na+ charges shown in Fig. [T] are quite 
distinct for the two Col sites (or the two Co2 sites), the 
electronic and magnetic properties only marginally dif- 
fer within each group of two sites. This establishes that, 
when dealing with the electronic properties, one might 
only consider a two site structure, with Co2 on a kagome 
planar structure, and Col on the complementary trian- 
gular structure. 



We have reached now a rather complete experimen- 
tal determination of the NMR parameters for the four 
Co sites formerly detected by NQR)^ which, together 
with the ^'^Na NMR spectra allowed us to determine the 
atomic structure of this Na2/3Co02 phascii^ The analy- 
sis of the NMR shifts and of the magnitude of the EFGs 
will allow us to confirm here the importance of cobalt 
charge disproportionation between the non-magnetic and 
magnetic sites, that is of the charge order induced in the 
structure. Wc shall discuss hereafter how this electronic 
kagome differentiation is reflected on the orbitals involved 
at the Fermi level at the different sites of the Co plane. 
The anisotropy of the EFG and NMR spin shifts will be 
seen to be related with the hole orbital ordering which 
governs the electronic properties of this phase. 



A. Site differentiation 

The data of Table [J allowed us to evidence that for 
the Cola and Colb sites the symmetry of the magnetic 
shift tensor K is axial {Kx=Ky), as found before for 
the EFG tensor (r; = or lyx—i^v)- The local symmetry 
of these sites had been a key argument used to resolve 
the crystal structureJ^ We also fully confirm here that 
the electronic properties are quasi identical for Cola and 
Colb sites which are slow relaxing and weakly sensitive 
to the magnetism. 

The knowledge of this actual atomic structure has al- 
lowed us also to resolve fully the difficulties connected 



B. Charge ordering 

This differentiation of sites had clearly led us to sug- 
gest that the ionic characters of the Col and Co2 sites 
arc not the samCfi^ We have suggested for long that 
the Col sites had an isotropic NMR orbital shift which 
is a strong signature that the Col sites are themselves 
nearly Co^^^, the orbital shift magnitude being further- 
more nearly identical to that found in the band insula- 
tor phase NaiCo02r^'^ in which all Co sites are non- 
magnetic Co'^"'". More accurately, as we are dealing with 
a metallic system the filling of the bands involving the 
Co t2g orbitals at the Fermi level is quite different on 
the Col and Co2 sites, and we can expect that the hole 
filling e of the t2g orbitals is small on the Col site. Let 
us assume that the hole content on Co2 sites is 6. As 
the number of holes per C0O2 formula unit is 1/3, the 
charge neutrality for the kagome structure implies then 
{3e + 9S)/12= 1/3, that is 

(5 == 4/9 - e/3. 



1. (a) Isotropic spin susceptibilities and Ti data 

The T dependent NMR shifts of the Co sites, which 
scale with that of ^'^Na give a measure of the spin sus- 
ceptibilities on the Co sites. If we consider in Table [IT] 
the isotropic shifts Kiso = (1/3) J^^aa on the Co sites. 
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we immediately find that the average values for the two 
Col sites Kiso{Gol)/^^K « 1.3 is an order of magnitude 
smaller than that for the Co2 sites Kiso{Co2)/'^^K w 
10.6, which indicates that the unpaired spins are domi- 
nantly on the Co2 sites. Similarly, both from NMR and 
low T NQR data we did find a ratio «70 between the spin 
contributions to 1/Ti for the two sites. Assuming similar 
hyperfine couplings for the two sites, one could antici- 
pate from these results a ratio of -s/TO « 8 of the local 
spin susceptibilities on the two sites, that is e/S ^ 1/8. 
This corresponds to an upper limit of 0. 06 for epsilon. 



where parameters 700 and Rei are the Sternhcimer anti- 
shielding factors, which characterize the enhancement by 
the core electrons of the atom of the EFG on the bare nu- 
cleus with respect to the EFG due to the outer electron 
distributions. 

The first contribution arises from all ion charges out- 
side the ion under consideration and can be calculated in 
an approximation assuming point charges 



K 



latt 



E q,{3cos^6i ~ 1) 
y,3 ' 



(11) 



2. (b) Isotropic orbital susceptibilities and orbital 
contribution to (TiT)^ 

Let us consider the difference found for the orbital 
shifts between the Col and Co2 sites. While for Col 
sites K"'^'' was found isotropic already in Ref. [ij, we do 
confirm here that K°^'' is anisotropic in plane for the 
Co2 sites with a sizable increase of K°^^. We can there- 
fore consider that the contribution of the hole orbitals to 
the orbital susceptibility is responsible for the deviation 
of K°^^ with respect to 1.91(3)% found in NaiCoOa-^ 
In a rough order of magnitude comparison we might con- 
sider then that the hole density on the Co sites is propor- 
tional to K°^^ — K°^^[x = 1), as we already suggested 
in Ref. [Tj. We do find then, from the results of Ta- 
ble |n]that i^°rj'(Col) - K°i;^{x = 1) = -0.01 ± 0.05(%), 
while K°rJ'(Co2) - K°^!;{x = 1) = 0.54(%). This would 
give again an upper estimate of e/5 < 1/10. Similarly 
knowing^^ that for Co'^"*' in insulating NaiCo02 the spin 
lattice relaxation rate is very weak, one could anticipate 
that the metallic orbital contributions to {TTi)~^ would 
be linked to the hole doping on the corresponding site. 
So there is no surprise to find out that {XT"'^'^)'^ is much 
larger for Co2 than for Col. 

So the data lead us to consider that we may neglect e 
as a first approximation and keep 5 ~ 0.44, although only 
accurate band calculations could tell exactly whether the 
partly filled bands at the Fermi level do not include any 
partial Col hole orbital in the kagome structure. 



3. (c) EFG induced by charge order 

Let us consider now the quadrupole frequencies vq^ 
which are directly proportional to the EFG on the nu- 
clear site probe. The latter arises from a non-symmetric 
distribution of electric charges around it. These charges 
can originate from non-bonding electrons, electrons in 
the bonds and charges of neighboring atoms or ions. The 
components of the EFG tensor can be written as the sum 
of two terms, the lattice Vj^^** and electron Vj^^ contribu- 
tions: 



K.a = (1 - 7oo)Kila" + (1 - Rel)V, 



*■ oca 7 



(10) 



where qi is the charge of ion i, located at distance tj of 
the probe nucleus, in a direction at an angle 9i from the 
Z axis of the EFG tensor. 

The second term in Eq.[TU] arises from unfilled electron 
shells of the orbitals of the considered site and of the 
distortions of the inner electron orbitals that they induce, 
embedded in the known value of (1 — Rei)- 

We have then to understand the relative weight of these 
two contributions in the measured EFG in this particular 
case of cobaltates. For a fully filled i2g multiplet such as 
the Co^"*" low spin state, the ionic structure is isotropic 
and the second term in Eq. [TU] should contribute negligi- 
bly to the on site EFG. 

Let us then point out that we can deduce significant 
indications about the respective weights of these two 
terms b y c onsidering the point charge calculations done 
in Ref. [2ll There we deduced average values for the 
Cola and Colb sites of vq = 1.32 MHz and a similar 
value VQ ~ 1.48 MHz for the average of the Co2a and 
Co2b sites. Experimentally the corresponding respective 
values VQ « 1.29 MHz and vq w2.36 MHz differ signif- 
icantly. This leads us to suggest then that the excess 
contribution of about 0.9 MHz on the Co2 sites is due to 
the on site hole orbitals. Furthermore, although we can- 
not rely fully on the numerical values of the point charge 
calculation, the agreement found with the Col data is 
indeed quite compatible with a negligible contribution of 
the on site orbitals on the Col sites, that is with a very 
small e value. 

To check the consistency of this analysis we may com- 
pare the data with that obtained for x = 0.35, for which 
we found that the Co sites are uniformly charged.— In 
that case, assuming a simple ordered distribution of Na2 
sites and x = 1/3, the point charge calculations give val- 
ues of Vq k, 2.5 MHz for all Co sites, whatever their 
location with respect to the Na sites, while the NMR 
data of Ref. [ij corresponds to vq ~ 4.1 MHz. So this 
much larger value would indicate that a 1.6 MHz on site 
contribution should be associated to 2/3 hole per Co site. 
Such a contribution to vq of 2.4 MHz/hole would corre- 
spond to 1.1 MHz for 0.44 holes on the Co2 sites in the 
X ~ 2/3 phase, remarkably close to the 0.9 MHz value 
estimated above. 

So, although these point charge estimates are certainly 
not fully reliable, they indicate that the excess EFG found 
on the Co2 sites has the right order of magnitude and is 
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compatible with the poor hole occupancy of the Col sites, 
as obtained here-above from the analyzes of K^''° and 



C. In plane anisotropics and orbital order 

The simple consideration of the atomic kagome struc- 
ture displayed in Fig. [11] allows one to expect an asymme- 
try of the local properties at the Co2 site, and to locate 
the two orthogonal axes {X*,Y*) which should be the 
Co2 site in plane local principal axes. 

Concerning the electronic structure, it is quite well 
known that the metallic bands of cobaltates are built 
from the Co t2g orbitals which subdivide in the aig or- 
bital which is axial and perpendicular to the plane and 
two subsets of in plane e' orbitals. From LDA calcu- 
lations of the band structure, done so far for uniformly 
charged cobaltate planes (see for instance Ref. |3^ , one 
usually expected a Fermi surface containing an aig sheet 
at the Fermi level centered at the F point and e' pockets 
at the zone cornersi^ Only the larger aig electron pocket 
has been observed so far by ARPES experiments, and 
that whatever the hole content j '^^'^° It has often been 
considered that, even in such simple uniform representa- 
tions of the C0O2 plane, the electronic correlations could 
largely influence the electronic structure by narrowing 
the bands and possibly pushing the e' hole orbitals be- 
low the Fermi level.— 

All that might not be relevant to the Na ordered and 
charge disproportionated structures, which would give 
a smaller Brillouin zone and a larger number of possi- 
ble bands at the Fermi level. But surprisingly here we 
find no experimental sign for multiband behaviour. In 
any case the dominant magnetic Curie- Weiss behaviour 
is governed by a single band in which correlations are im- 
portant. So far the nature of the Co orbitals which par- 
ticipate in this correlated electronic band are not readily 
available from the data. It remains then important to 
try to determine experimentally what are the respective 
weights of the aig and e' orbitals involved at the Fermi 
level in the actual electronic structure of the 2/3 phase. 

We consider that the non axial NMR parameters deter- 
mined here may help to answer this question. Assuming 
that Col is in a completely filled shell Co^+ ionic state 
(that is e=0), if the electronic structure only results from 
aig hole orbitals on the Co2 sites, the non axial NMR pa- 
rameters would only be induced by the absence of trans- 
fer integral between Co2 and Col sites. But if a fraction 
of the hole density on Co2 resides on e' orbitals, the local 
asymmetry of the transfer integrals would induce distinct 
occupancies of the two e' orbitals on the Co2 sites. One 
could expect a larger hole occupancy of the e' orbital 
pointing in the X* direction in Fig. 1111 which would cor- 
respond to an ordering of the in plane hole orbitals on 
the Co2 sites. In any case such an unbalanced popula- 
tion of the e' orbitals would induce an on site anisotropy 
of the NMR parameters on Co2 sites. We shall here try 




FIG. 11. (Color online) The kagome structure in the C0O2 
planes allows to clearly locate two orthogonal axes {X* , Y") 
which should be the Co2 site principal axes in the plane. 



to estimate the aig and e' hole populations by consider- 
ing first the asymmetry of the EFG and then that of the 
NMR shifts. 



1. In plane EFG anisotropy on Co2 sites 

Both asymmetries of the distributions of Na ions and 
Co charges located in the plane do introduce an in plane 
EFG asymmetry on the Co2 sites. We could evaluate 
the part of this asymmetry due to ionic charges from the 
point charge calculations of the first term of Eq. [TU] done 
in Ref. [2l|. We indeed found that the distribution of Na+ 
charges induces a slight tilt of the Z principal axis of the 
EFG with respect to the c axis on both Co2 sites. The 
overall ionic charge distribution also results in X and Y 
principal axes of the EFG which are distinct on the Co2a 
and Co2b sites and which were found markedly modified 
if we displace the Na atoms with respect to their ideal po- 
sitions as suggested by the x-ray Rietveldt data analysis 
of Ref. [2l|. However in all cases considered, the computed 
values of | Vyy — Vxx \ was found to correspond to a dif- 
ference of NQR frequencies (^'y — vx) hi^ a range between 
0.4 and 0.9 MHz. This appears on the average signifi- 
cantly smaller than the 0.8 and 0.9 MHz values deduced 
from the experimental asymmetry parameter 77 sa 0.35 
found for both Co2a and Co2b sites (see Table H]). 

We may then suggest that some of the non axial charge 
distribution is linked with the on site term. As the aig 
orbital is axial, any on site non axial contribution to the 
EFG of the Co2 could only result from an unbalance in 
the populations of the e' orbitals4^i^. The fact that the 
asymmetry parameter and EFG on the two Co2 sites are 
quite similar further supports the existence of a sizable 
on site e' hole orbital contribution to the EFG. Unfor- 
tunately the large sensitivity of the point charge calcu- 
lations of the EFG to the actual atomic structure does 
not allow us to subtract reliably this EFG contribution 
from the experimental data, and to deduce then the prin- 
cipal axes and magnitude of the on site contribution to 
the EFG. It is clear then that more accurate EFG deter- 
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minations based on ab initio band calculations are then 
needed to relate the EFG data to the relative populations 
of the flig and e' orbitals. 



2. In plane Knight shift anisotropics 

While the anisotropy of Na charge distribution con- 
tributes to the EFG anisotropy, the Na ions are not in- 
volved directly in the magnetic properties, which arc fully 
associated with the C0O2 layers. So, the anisotropies of 
the NMR shifts of the Co2 sites should give much better 
indications on the anisotropy of the in plane electronic 
structure. The non axial values of the Co2 shifts summa- 
rized in Table Hi] are indeed pointing out the existence of 
an in plane electronic anisotropy. One can notice there 
that both 7^°''^ and K^P'^'iT) are much larger in the X 
direction than in the Y direction. Let us note that these 
experimental directions are those identified from the EFG 
tensor, with the usual convention \Vyy\ > I'^xxl, so that 
the smallest quadrupole frequency occurs in the planar 
direction for which the NMR shifts are largest, as seen 
directly in the spectra. However so far we do not know 
for sure how {X,Y) relate to the {X*,Y*) axes defined 
in Fig. [m 

We did find as well an anisotropy of the spin contri- 
butions to (TTi)~^ which correlates with the anisotropy 
of hyperfine couplings, so that no really new informa- 
tion is embedded in this result, especially on the in plane 
anisotropy. 

a. Orbital shift in plane anisotropy SQUID data 
had allowed us to evidence that the macroscopic suscep- 
tibility is larger along the ab plane than in the c direction. 
This data allowed us to establish that the T dependence 
of the susceptibility, which is determined by the spin con- 
tribution ^'^P'" (T) is nearly isotropic, which has been 
confirmed as well from ^"^Na NMR data^^i So the sus- 
ceptibility anisotropy is dominated by that of x°^'' and 
Xab'' > Xc^'' is sufficiently large to be practical to align 
the single crystal grains in the applied field at room tem- 
perature. But these SQUID data on powder crystallite 
samples do not give any indication about the in plane 
anisotropy of orbital susceptibility. The present NMR 
data in Table |TT] establishes a planar anisotropy of x°^^ 
on the Co2 sites, for which K'^'' > K^'''' > K^''''. 

Let us point out that K°^^ can be a purely ionic term, 
as is the case for Co'^"'", but states at the Fermi level might 
also give a specific contribution. This is indeed expected 
as we found such a metallic band orbital contribution 
to (TiT)-i for both Col and Co2 sites. In fact Xf '' 
is unmodified with respect to that on Co'^"'", so that the 
extra contributions found for K'^'' and Ky'' could be 
mainly due to such band contributions, as is suggested 
by the correlation between the spin and orbital in plane 
anisotropies found in Table HIl One could notice however 
hereagain that if the hole orbitals were purely aig, one 
would not expect any in plane anisotropy of K°^'', so 
that heregain band calculations of x^a niight help to 



determine the populations of the aig and e' orbitals. 

These observations on the Na and Co sites are remi- 
niscent of the respective situations found for the NMR 
of the ^^Y and of the ^'^Cu nuclei of the Cu02 planes of 
the YBCO cuprates. In that case, as for Na NMR here, 
the ^^Y NMR shift which is only coupled to spin mag- 
netism was found nearly isotropics^ while the anisotropy 
of orbital shift of *''^Cu was due to that of x°^''i though 
in that case the orbital contribution due to the d^2_y2 in 
plane hole orbital^ does not display any anisotropy in 
the Cu02 plane. 

b. Spin shift in plane anisotropy As the spin sus- 
ceptibility is nearly isotropic from SQUID and Na NMR 
shift data, the large in plane anisotropy of the Co2 NMR 
spin shift data K^p™{T) has to be associated with an 
anisotropy of hyperfine couplings, which are linked with 
the occupancy of the electronic orbitals involved at the 
Fermi level. 

Further comparison with the cuprates could help us to 
better understand the in plane anisotropy of ^^Co NMR 
spin shifts detected here. Quite generally the spin shift 
tensor K^J^''^ is slightly more complicated than indicated 
in Eq. [5] as nuclear spins at site i can be coupled both 
to the on site magnetism and to that of their j near 
neighbour sites, so that 



^aa (v — ^a, 
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where Aaa is the anisotropic on-site hyperfine coupling 
tensor and Baa is the transferred hyperfine coupling 
through oxygen orbitals with the j nearest neighbors. In 
YBCO there is no hole occupancy on the Y site and there- 
fore no on site hyperfine coupling Aaa so that the hyper- 
fine coupling of ^^ Y to the Cu site magnetism is only due 
to the transferred Baa which has been found isotropic 4^ 
On the contrary the anisotropy of on site coupling Aaa 
is responsible for the large measured (c, ab) anisotropy 
of the ^'^Cu spin shift. Even in this simple case with a 
single Cu axial site it took sometime to understand that 
a rather large isotropic value of Baa was required to ex- 
plain the data4^ Here, we shall benefit of this approach 
to perform a tentative preliminary analysis of the Col 
and Co2 data summarized in Table HIl 

For the axial Col sites, only a (c, ab) shift anisotropy is 
detected. If e=0, then Col sites are non-magnetic Co'^"'" 
ions which bear no on site magnetic spin susceptibility. 
As the Y ions in the cuprates they would be inert magnet- 
ically and would sense the magnetism of the Co2 orbitals 
only through transferred hyperfine couplings Baa with 
the Co2 magnetism. The (c, ab) anisotropy should then 
be associated with an anisotropy of the coupling Baa of 
the Col with its first Co2 nearest neighbours. A dipolar 
coupling of Col with its six Co2 neighbours would cor- 
respond to a larger negative contribution to the shift for 
H II c, contrary to observation. A non zero e value with 
a small filling of the aig hole orbital on the Col would 
better explain the sign of the (c, ab) anisotropy with a 
maximum shift contribution for -ff || c. 
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But of course the most intriguing experimental infor- 
mation resides in the totaUy non-axial spin shift values 
found here for the Co2 sites, given in Table |lll Aver- 
aging the data for the Co2a and Co2b sites results in 
Kf"'/^^K = 18.9; K^P'^'/^^K = 7.5 and i^f ™/ ^^K = 
4.9. 

These large anisotropic contributions to the shift might 
be dominated by on site hyperfine couplings with the aig 
and e'g orbitals. We can note hereagin, as for the discus- 
sion of the EFG, that the axial symmetry of the aig which 
points along the c axis would only yield axial shift contri- 
butions and would not differentiate X and Y directions, 
contrary to the e'g orbitals. So the large difference be- 
tween Kx and Ky can only be attributed to the in plane 
order of e' orbitals associated with the kagome structure. 
If one considers the transfer integral paths of e' orbitals 
between Co2 neighbouring sites, one would be inclined 
to associate the X direction, with the largest shift value 
with the direction X* sketched in Fig. [TTI 

Let us note however that nothing forbids as well to 
consider here that the Baa term itself is anisotropic. De- 
tailed consideration of the Wannier orbitals constructed 
for the full atomic structure of this ordered x = 1/3 phase 
could help to explain these measured anisotropics and 
would allow altogether to determine the relative popula- 
tions of aig and e' hole orbitals on the Co2 sites at the 
Fermi level. 



VII. CONCLUSION 

In this experiment we have completely determined ex- 
perimentally the NMR parameters of the four Co sites 
of the atomic structure of this Na2/3Co02 phase, which 
allows us to confirm that its electronic structure differ- 
entiates in fact two sites in the C0O2 plane: Col with a 
charge 3 + e and Co2 with a charge 3.44 — e. The data 
for the EFG, the orbital NMR shifts and orbital TiT are 
compatible with a hole content at most of e < 0.06 on 
the Col site as we had already anticipated. 

From our NMR data we confirm as well that a single 
band dominates the anomalous electronic properties of 
this phase at the Fermi level. A totally unexpected aspect 
revealed by these data is the existence of large in plane 
anisotropics of the EFG and NMR shifts on the Co2 site. 
We have shown that this implies that the hole orbitals 
on those sites are not exclusively on aig orbitals and that 
parts of the holes reside in the e/g in plane orbitals which 
should then play an important role in the delocalization 
of the holes. 

So the single hole band hybridizes those orbitals and 
the kagome like electronic structure that we evidenced 
then displays not only a charge disproportionation be- 
tween Col and Co2 sites, but as well an in plane charge 
order of the t2g orbitals on the Co2 sites. Quantum chem- 
istry calculations taking into account the full structure 
might permit to go beyond the qualitative arguments de- 
veloped here and should provide quantitative data for the 



t2g energy levels from the values obtained experimentally 
for the EFG, hyperfine couplings and orbital shift param- 
eters. 

One would like then to understand the respective roles 
of the Na order and of the electronic correlations in driv- 
ing this electronic structure. In principle LDA calcula- 
tions including the Na order should allow to understand 
whether the Na potential is sufficient to shift down the 
Col t2g energy levels and to induce the Co^^^ filling of 
the Col orbitals. Preliminary calculations along this line 
seem to indicate that this is not the ease and therefore 
that it is essential to better take into account the elec- 
tronic interactions 41 In such a case the kagome like orga- 
nization of the electronic structure would be an intrinsic 
disproportionation in the Co planes which would lock the 
Na order 4^ 

Our work evidences that the resulting situation is quite 
simple in this Na2/3Co02 phase, and resumes in a sim- 
ple two site charge disproportionated C0O2 plane which 
could be introduced in a simplified model. This would 
allow one to perform then cluster DMFT calculations to 
establish the role of correlations which are of course es- 
sential to explain the large Curie- Weiss like T dependence 
found for the susceptibility as well as the 2D ferromag- 
netic correlations deduced from our spin lattice relax- 
ation dataJ^ 

In any case we have shown that the Nai:Co02 com- 
pounds display stable Na ordered phases in which charge 
ordered states of the C0O2 planes occur so that the real 
atomic structure has much larger 2D unit cells than that 
of the C0O2 plane. We might anticipate then that most 
ARPES experiments, which isolate a specific surface after 
cleaving do not sample the electronic properties of the ac- 
tual 3D ordered state but are affected by the unavoidable 
Na disorder in such experimental conditions, although in 
some cases specific local Na orderings have been observed 
by STM experiments on cleaved surfaces4^ 

An interesting case is that of some misfit cobaltates^ 
for which the multilayer structure is obtained by stack- 
ing rocksalt structures with C0O2 planes. In that case 
cleaving is done in the middle of the rocksalt structure, 
as in High Tc cuprates, and an ARPES experiment^ has 
evidenced that in the specific case of BaBiCo02 the elec- 
tronic structure displays a reconstruction which could re- 
sult from a charge order analogous to that seen in the 
Na^Co02 phases. Unfortunately in that case it has so far 
not been possible to perform NMR experiments allowing 
to resolve the sites?^ so that detailed comparisons are 
not possible. 

Finally we have established here on a specific case the 
importance of charge disproportionation and orbital or- 
dering in this compound. No doubt that for all hole 
dopings for a; > 0.5 similar effects do occur as diverse 
Na orderings are anticipated and observed and result 
in very distinct correlated electronic states . ^'""^"'^^'^^ We 
think that the present NMR studies allow us to reach for 
the first time a very detailed description of such effects 
in multiorbital transition oxides. 
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Appendix A: Oxygen vacancies 

A recent publication^^ is dedicated to discredit our 
resultsi^i'^ and in support of structural models of the Na 
cobaltates that this group has been proposing for many 
years for large x values. 

(i) In this paper the authors suggest that our sam- 
ples, being synthesized at high T, are not Na2/3Co02 but 
rather Na2/3CoOi.98. We can hardly discard solely from 
our NMR data the presence of such a low content of oxy- 
gen vacancies. As each Co has six near neighbor oxygen, 
only about 6% of the Co would have such near neigh- 
bor vacancies. If those sites would get well defined NMR 
shift and EFG they should be detected either in NMR or 
NQR. But as soon as some disorder is expected in the cor- 
responding parameters due to long distance interactions 
between defects, we would not be able to detect them. 
We may however easily exclude that vacancies could in- 
duce perturbations beyond their first nearest neighbors. 
Those would yield a broad background signal involving 
then a large Co NMR intensity. So, if present, such oxy- 
gen vacancies should just play a similar role to that of Na 
atoms and give delocalized carriers on the kagome struc- 
ture, independently of the Col sites Co'^^ signal which 
are connected with the ordered Na sites. 

To conclude on that point, there is so far no mi- 
croscopic experimental evidence for the existence of 
such oxygen vacancies neither in Ref. |2^ nor in our 
NMR/NQR data. 

(ii) Beyond this materials research aspect, Shu et 
al. speculate, heragain without any microscopic ev- 
idence, that the Na order in our samples (puta- 
tively Na2/3CoOi.98) might be a structure built from 
a succession of planes with Na trivacancies and Na 



tetravacancies^ which they have proposed for their sam- 
ple with an assumed composition of Nao.7iCo02. 

We claim that this is definitely contradicted by all 
our former NMR, NQR and x-ray dat a^^i^^ and by the 
present ^^Co NMR data. It is easy to establish that this 
structure corresponds to 5 Na crystallographic sites (with 
relative intensities 3/6/3/2/3), and 8 Co sites (with rel- 
ative intensities 1/1/1/3/3/3/3/9). This is well beyond 
the three Na sites we have detected by NMR^ and the 
four Co sites we have detected both by NQRiSii^ and 
here by NMR. Finally the structure we have proposed 
from the analysis of our NMR/NQR data in Ref. [l^ has 
been confirmed by the perfect fit obtained with the Ri- 
etveldt x-ray data analysis performed on our samples<^ 

So the self consistency of our data on our Na2/3Co02 
samples cannot be put into question. Whether the struc- 
ture proposed by Shu et al. for their Nao.7iCo02 sam- 
ples is the actual one remains a question which might 
be resolved in the future by performing independent 
NMR/NQR investigations on their samples. 



Appendix B: Check of the pow^der alignment by 
NMR 



The alignment of the sodium cobaltates powders being 
dependent on the microstructure of the grains, some of 
our samples were indeed not well oriented. We demon- 
strate here how we could use the NMR spectra to distin- 
guish poorly oriented samples from well oriented ones. 

In Fig. [12] we show a ^'^Na NMR spectra measured in 
two samples of Na2/3Co02 compound. One of them is 
well oriented sample - Sample A - and we show there the 
spectra measured in two directions of the applied mag- 
netic field - Ho -L c and Hq \\ c. As evidenced beforo^i^ in 
the Na2/3Co02 phase there are three unequvalcnt sodium 
sites - Nal, Na2a, Na2b. We found that although the 
central lines and satellites singularities of the ^■^Na spec- 
trum can be easily located, the overall shape of the NMR 
spectrum is very sensitive to the degree of orientation of 
the powder sample. In Fig. [T^lwe also show results of the 
computer simulations of the ^'^Na NMR spectra using the 
parameters of Eq. [2]from Ref. 19). Here we have optimized 
the particles c axis angular orientation distribution to 
fit accurately the spectrum (Lorentzian distribution of 
±5° width). It's important to mention that we used ex- 
actly the same angular distribution parameters for the 
^^Co simulations shown in Figs. [31111 As one can sec the 
agreement between experimental and simulated spectra 
are very good. 

At the bottom of Fig. [H we show a ^^Na NMR spec- 
trum measured in sample B of Na2/3Co02, which hap- 
pened to be less well oriented. Comparing this spectrum 
with that measured in Sample A it is easy to see for 
Ho II c, the presence of Nal singularities specific to the 
Ho -L c direction of applied magnetic field. That is due to 
a fraction of particles which remained randomly oriented 
as they were either not single crystallites, or were too 
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FIG. 12. (Color online) ^^Na NMR spectra: comparison of 
oriented sample A (black open circles for Ho J- c and open 
squares for Ho \\ c) and sample B (red closed squares for 
Ho |[ c). The solid colour lines display the computed spectra 
for the three different Na sites, and the total spectum which 
fits perfectly the data for sample A. The presence of the Nal 
singularities specific for the Ho -L c direction spectrum of 
sample A in the Ho \\ c spectrum of sample B is clear and 
allowed us to conclude that sample A is better oriented than 
sample B. The extra line near the central transition at about 
5.88 T with small Knight shift has a very weak intensity, 
which corresponds at most to 5% of all Na nuclei, and might 
be assigned to defect sites in the ordered Na structure, or to 
a minute amount of impurity phases present in these samples 
with large crystallites. 



much packed or agglomerated and could not orient. Also 
a less prominent, but still detectable feature of this less 
oriented sample is the poor resolution of the quadrupolar 
satellites. 

In Fig. [13] we demonstrate the distinction between the 
^^Co NMR spectra measured in these two samples as 
well. The comparison can summarized as follows 

1) First of all in the sample B spectrum measured for 
Ho II c one can see an extra line, circled as (1) which 
corresponds to the central lines of Cola and Colb sites 
in the Hq _L c direction - marked by arrow in Fig. 1131 

2) Also one can see that for the well oriented sample 
A that the NMR spectrum measured for Hq || c is rather 
symmetric whereas the spectrum for Hq _L c is asym- 
metric with a higher signal intensity for the lower field 
values which could be assigned to fast relaxing cobalts. 
In sample B even in the Hq \\ c spectrum some asym- 
metry could be seen and is here again associated with 
the contribution from grains which did not orient in the 
field. This extra intensity, circled as (2), had mislead us 
by suggesting that a Co site labeled as Co3 could have a 
large purely axial shift (see Fig. 2c in RefJ^). 

3) Again in the sample B spectrum the quadrupolar 
satellites are less resolved as one can see in Fig. [T31 on 
the part circled as (3). 



FIG. 13. (Color online) Comparison of ^®Co NMR spectra 
of the well oriented sample A (black open circles for Ho -L 
c and open squares for Ho || c) with that of the less well 
oriented sample B (red closed squares for Ho \\ c). The extra 
contributions in the Ho \\ c spectrum of sample B due to the 
Col and Co2 signals of non-oriented particles are highlighted 
within circles (1) and (2). 



Appendix C: Determination of K{T) from full 
analyzes of the ''^Co NMR spectra 

The ^^Co NMR of the slow relaxing Col sites in the 
Na~o.7Co02 compounds are easy to observe and were 
seen even by continuous wave NMR.— In the Na2/3Co02 
compound even small details for the Cola and Colb sites 
could be detected and have been reportedi^ The axial 
symmetry of these sites symmetrizes the magnetic shift 
and quadrupolar tensors which makes the NMR lines of 
^^Co in these positions rather narrow and intense. 

The case is more complicated for the fast relaxing Co2 
sites in the Nafso.7Co02 compounds. Even in the single 
crystals the NMR spectra of these sites are broad and 
not well resolved 1^ The same situation occurs in our ori- 
ented powder samples - see Fig. [5] as an example. As in 
our samples the c-axis of different crystallites are aligned 
along Hq II c one can distinguish in the spectrum the two 
sets of 7 lines which correspond to the different transi- 
tions TOO- (m — 1) for the two cobalt sites Co2a and 
Co2b. For Hq _L c, the NMR spectrum of the fast re- 
laxing cobalts is found broad and non-resolved, as the 
particles in our samples have a powder distribution of 
ah axes for Hq _L c. This experimental fact immediately 
tells that the local symmetry for the Co2 sites is lower 
than for the Col sites. 

To clarify the spectra of the Co2 sites for Hq _L c we 
used various facts: 

1) The position of the different singularities in the fast 
relaxing Co2 spectra are the same whatever the sample 
and whatever the sample preparation method, as one can 
see in the Fig. [21 The positions of the singularities in the 
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NMR spectra depend on both the quadrupole sphtting 
(which could be considered practically as temperature 
independent) and the magnetic shift which scales with 
the susceptibility Xs(T) . 

2) The main feature of the Na2/3Co02 phase is a large 
temperature dependence of the spin susceptibility. As 
we have shown that the local fields are quite similar for 
the three detected Na sites^, the first moment (center of 
gravity) of the ^^Na. NMR spectrum 23^iso ^ Axs{T) 
allows us to follow the spin susceptibility Xs{T) of the 
cobalts up to room T. Therefore using ^^Kl^°{T) as a 
reference we tried to follow the corresponding T variation 
of the Co shift tensor. 

We measured then the ^^Co NMR spectra at different 
temperatures in swept field mode and at each tempera- 
ture pointed the positions (field values) of all singularities 
in the spectra of fast relaxing Co2. After that we con- 
verted them into field shift values using Eq.Uand plotted 
them versus '^^K{T) taken at the measurement temper- 
ature. Each spectrum taken at a given T then allows to 
report all the data points on a vertical line in Fig. [Tik . 

As one can see, some singularities for the fast relaxing 
Co2 appear well organized around extremal values of the 
shift, with 3-|-3 satellites spaced by the quadrupole split- 
ting Va - in Fig. 114b we show two such sets of singularities 
by lines. Initially we assumed that these two sets of sin- 
gularities correspond to the two cobalt sites: one with 
a large shift in Hq ± c direction and the other a small 
shift, but this assumption did not allow us to explain all 
the data. 

Unexpectedly we found that the values of quadrupole 
splitting Va for these two sets of singularities correspond 
to the i^x and i/y values of the Co2b site (the correspond- 
ing i/Q and ?7 values were determined by NQR with high 
presicion^i) . Therefore the main conclusion of Fig. [Hb 
is that we have a large anisotropy of in-plane magnetic 
shift for the dominant Co2b site with distinct values of 
Kxx and Kyy- 

As the number of Co2a sites in the unit cell of the 
Na2/3Co02 phase is twice smaller than that of Co2b sites, 
the intensity of their NMR signal is also twice smaller. 
This makes more difficult to resolve the Co2a site in the 
Hq _L c spectra. As one can see in Fig. [Hb many sin- 
gularities are located in between the identified Co2b sin- 
gularities. Using them and the values of vq and rj for 
Co2a sites)^ one can suggest two sets of singularities for 
the Co2a sites - see Fig. [Hb . And in fact the Co2a sites 
exhibit the same in-plane magnetic shift anisotropy as 
the Co2b sites - just with slightly smaller values of the 
magnetic shifts and quadrupole splittings. 

Therefore the fast relaxing cobalt Hq _L c spectra can 
be finally explained in a simple manner. In Fig llSI we 
show the simulated set of NMR lines of Co2b for the two 
outermost cases for which Hq \\ X and Hq \\ Y. These 
two spectra form the main singularities in the cobalt 
spectra while most of its intensity comes from the parti- 
cles with intermediate orientations. The large difference 
in the Kxx and Kyy values makes the cobalt Hq _L c 
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FIG. 14. (Color online) (a) The global summary of all our 
magnetic shift data for the various singularities of the fast 
relaxing cobalt spectra, plotted versus the isotropic compo- 
nent of the ^'^Na NMR shift. The most pronounced singulari- 
ties in the experimental spectra are marked by filled symbols, 
whereas weaker singularities are displayed with open sym- 
bols. DilTerent symbols correspond to spectra for differently 
prepared samples of the same phase, (b) Here, the same data 
are reported, but with lines drawn which underline the central 
line and 6 quadrupolar satellites of ^®Co in the Co2b positions 
for Ho II X (blue full lines) and Hq \\ Y (green dashed lines) 
directions; (c) On the same data again, we have drawn the 
seven ^^Co NMR lines which correspond to the Co2a posi- 
tions for Ho II X (red full lines) and Hq \\ Y (magenta dashed 
lines) directions. 
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FIG. 15. (Color online) The ^^Co NMR spectrum of fast re- 
laxing Co2 is shown by full dots joined by a continuous line. 
The positions of the seven lines estimated from the full data 
analysis are shown for Ho \\ X (blue solid line) and Ho \\ Y 
(green dashed line). The correspondence with the singular- 
ities in the experimental spectrum is clear. The position of 
these two set of singularities is then determined by these ex- 
tremal lines while most of the intensity in the spectrum comes 
from the particles with intermediate orientations. The large 
singularities around 6.4 Tesla comes from imperfect compen- 
sation of the large signal of the slow realaxing Col sites. 
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